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We show that the black hole laser effect discovered by Corley & Jacobson should be described in 
terms of frequency eigenmodes that are spatially bound. The spectrum contains a discrete and finite 
set of complex frequency modes which appear in pairs and which encode the laser effect. In addition, 
it contains real frequency modes that form a continuous set when space is infinite, and which are 
only elastically scattered, i.e., not subject to any Bogoliubov transformation. The quantization is 
straightforward, but the calculation of the asymptotic fluxes is rather involved. When the number 
of complex frequency modes is small, our expressions differ from those given earlier. In particular, 
when the region between the horizons shrinks, there is a minimal distance under which no complex 
frequency mode exists, and no radiation is emitted. Finally, we relate this effectto other dynamical 
instabilities found for rotating black holes and in electric fields, and we give the conditions to get 
this type of instability. 

PACS numbers: 04.70.Dy, 43.35.+d 



I. INTRODUCTION 



In 1998, Corley and Jacobson made the following interesting observation pQ. They noticed that the propagation 
of a superluminal dispersive field in a stationary geometry containing two horizons, as it is the case for nonextreme 
charged black holes, leads to a self amplified Hawking radiation (for bosonic fields). The origin of this laser effect 
can be attributed to the closed trajectories followed by the negative Killing frequency partners of Hawking quanta. 
Because of the superluminal character of the dispersion, the partners indeed bounce from one horizon to the other. 

Besides charged black holes, the acoustic geometries associated with a moving fluid that crosses twice the speed of 
sound also contain a pair of horizons. For one-dimensional flows these geometries are of the form 

ds 2 = -dt 2 + (dx - v(x)dt) 2 . (1) 

This expression obtains [2j|3] when considering the propagation of low frequency phonons in a moving fluid of velocity 
v(x). The sound velocity is assumed to be constant and has been set to 1. This metric possesses a black hole (BH) or 
a white hole (WH) horizon when the gradient d x v at the horizon \v\ = 1 is positive or negative. When the dispersion 
is superluminal (anomalous) all such systems should experience a lasing effect. This should apply in particular to the 
BH-WH pair realized in a Bose Einstein condensate [4]. 

The original analysis and that of [5 were both carried out using wave packets. In the present work we show 
that there exists a more fundamental description based on frequency eigenmodes which are asymptotically bound (in 
space). When v(x) is constant and subsonic for both x —> ±oo, there is a discrete set of complex frequency modes 
and a continuous family of real frequency modes. (If instead the subsonic region is finite and periodic conditions 
imposed, the situation is more complicated [7] and will not be considered here.) In our case, the real frequency 
modes are asymptotically oscillating and normalized by a delta of Dirac. They are not subject to any Bogoliubov 
transformation. In fact, the scattering matrix at fixed uj only contains reflection and transmission coefficients mixing 
right and left moving (positive norm) modes. This was not a priori expected since the matrix associated with a single 
BH (or a WH) is 3 x 3 and mixes positive and negative norm modes in a nontrivial way |8j |9] . 

The discrete set is composed of modes that vanish for x — > ±oo. They form two modes subsets of complex conjugated 
frequencies, each containing a growing and a decaying mode. The time dependence of the coefficients in each subset 
corresponds to that of a complex upside down and rotating harmonic oscillator. Both the real and the imaginary part 
of the frequency play important roles in determining the asymptotic properties of the fluxes. We notice such modes 
were encountered in several other situations [10-15 . We also mention that a stability analysis of BH-WH flows in 
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Figure 1: Velocity profile v(x) as a function of kx, for Kb = k>w, D — 0.33, and kL — 8. The horizons are located at kx — ±8, 
where v(x) + 1 = 0. 



Bose Einstein condensates was presented in [16]. We reach different conclusions because we use different boundary 
conditions. 

In Sec. [IlJ we present our settings. In Sec. we demonstrate that the set of spatially bound modes contains a 
continuous part and a discrete part composed of complex frequency modes. In Sec. |IV| and Sec. [V| we study the 
properties of the modes, and show how the complex frequency modes determines the fluxes. We also relate our 
approach to that based on wave packets [IJ[5], and explain why the predictions differ in general, and in particular 
when the number of discrete modes is small. In Sec. [VI] we give the conditions to get complex frequency modes in 
general terms. 



II. THE SETTINGS 



We work in two space-time dimensions and consider the stationary metrics of Eq. ([TJ which contain a BH-WH pair. 
We restrict ourselves to flows that are asymptotically constant, i.e. we consider velocity profiles such as 



v(x) 



-1 + D tanh 



n B (x - L) 



D 



tanh 



D 



(2) 



see Fig. [I] The BH (WH) horizon is situated at x = L (x = —L). We suppose that the inequality nL/D ^> 1 
is satisfied for both values of ft, where D g]0, 1]. In this case, the two near horizon regions of width Ax ~ D/k 
are well separated, and the surface gravities of the BH and the WH are, respectively, given by Kb = O x v\ x= l and 
Kw = —dx v \x=-L- The minimal speed \v- \ < 1 is reached for x — >• d=oo, whereas the maximal speed > 1 is found 
at x = between the two horizons. When nL/D ^> 1, their values are 



v± 



(3) 



As emphasized in [8, 9], when using nonlinear dispersion relations, D fixes the critical frequency cj m ax above which 
no radiation is emitted by a single BH, or WH. Similarly here, there will be no unstable mode above co> m ax- 
As in [1 we work with a real field <p obeying a quartic superluminal dispersion relation 



k 2 + k 4 /A 2 . 



(4) 



The UV scale characterizing the dispersion is given by A, and ft is the frequency measured in the preferred frame, 
that comoving with the fluid. Most of the results we shall derive also apply to higher order superluminal dispersion 
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relations. They also apply to the Bogoliubov theory of phonons in Bose condensates which is a variant of the present 
case, as can be verified by comparing [9] to [8]. The action of <\> in the metric of Eq. is 



dtdx 



(d t <p + vd x 4>) 2 -(d x <p) 2 -^(d 2 x 4>) 2 



and the wave equation is 



(d t + d x v)(d t +vd x )-d 2 x + ^dt 



= 0. 



(5) 



(6) 



When the flow is stationary, one can look for solutions with a fixed frequency A = idf. Inserting 
Eq. d6b yields 



i\t 



+ d x v)(-i\ + vd x ) - d 2 + 



0A=O. 



<j>x{x) in 



(7) 



Because of the quartic dispersion, the number of linearly independent solutions is four. It would have been n if the 
dispersion relation would have been Q 2 = p 2 + p n /A n ~ 2 rather than Eq. Q. However, when imposing that the modes 
be asymptotically bound for x — > ±oo, the dimensionality is reduced to 2 or 1 depending on whether A is real 
or complex, but irrespective of the value of the power n. (To avoid confusion about the real/complex character of A 
we shall write it as A = uo + iT, with u and T both real and positive. The other cases can be reached by complex 
conjugation and by multiplication by —1.) 

The necessity of considering only asymptotically bound modes (ABM) comes from the requirement that the ob- 
servables, such as the energy of Eq. (15), be well defined. Returning to Eq. (J5|, the conjugate momentum is 



d t <t> + vd x </>, 



(8) 



the scalar product is 



(01 102) = i dx [0i7T 2 - </>27Ti] 



and the Hamiltonian is given by 



H 



dx 



(d t <t>f + (i-v 2 )(d x( j>f 



A 2 



{d 2 x 4>f 



(9) 



(10) 



In the subspace of ABM, the Hamiltonian is hermitian, i.e. {(j)\\H(j)2) — (Hcfrilfa). 

We conclude with some remarks. First, when written in the form Eq . (|9| ), the scalar product is conserved in virtue 
of Hamilton's equations, and the hermiticity of H. Second, from Eq. (10), one sees that the Hamiltonian density is 
negative where the flow is supersonic: v 2 > 1. We shall later see that the supersonic region should be 'deep' enough 
so that it can sustain at least a bound mode, thereby engendering a laser effect. Third, when considering ABM, 
eigenmodes characterized by different frequencies are orthogonal in virtue of the identity [TO], [13] 



(A'-A*)(^v 



0, 



(11) 



which follows from the hermiticity of H. Finally, we notice that complex frequency ABM can exist in the present 
case because neither the scalar product Eq. ([9|, nor the Hamiltonian Eq. (10) is positive definite, see [10] p. 228. On 
the contrary, since fermionic fields are endowed with a positive scalar product [17 , no complex frequency ABM could 
possibly be found in their spectrum [T8] , 



III. THE SET OF ASYMPTOTICALLY BOUND MODES 



A. Main results 



In the BH-WH flows of Eq. (|2|, the set of ABM contains a continuous spectrum of dimensionality 2 labeled by 
a positive real frequency cj, and a discrete spectrum of TV < oo pairs of complex frequencies eigenmodes. We shall 
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suppose that this set is complete. 1 That is, any solution of Eq. Q with {<j)\4>) < oo can be decomposed as 

duj (e^ + a W) ,CW] + h.c.) 

(e" zAa£ b a <p a (x) + e~ iX ^ c a ^ a (x) + ft.c.) . (12) 

For flows that are asymptotically constant on both sides of the BH-WH pair, we shall show that the real frequency 
modes can be normalized according to 

= = o, (13) 

where the discrete index i takes two values u,v, and where uj^uj' > 0. The index u : v characterizes modes which are 
asymptotically left (v) or right moving modes (u) with respect to stationary frame. 

When A is complex, the situation is unusual. Yet it closely corresponds to that described in the Appendix of [TO] . 
In fact whenever a hermitian Hamiltonian possesses complex frequency ABM, one obtains a discrete set of two- modes 
(tpa^a) of complex conjugated frequency A a , A*. Their 'normalization' can be chosen to be 

(VV|V?a) =0, ftM^a) = iSa,a>: (14) 



with all the other (independent) products vanishing in virtue of Eq. (11). Since the overlap between modes belonging 
to the continuous and discrete sectors, such as {(j) u \^fa)^ a ls° vanish, eigenmodes of different frequency never mix. 
Moreover, since the positive norm modes (j) 1 ^ all have uj > 0, one cannot obtain Bogoliubov transformations as those 
characterizing the Hawking radiation associated with a single BH (or WH). This implies that the (late time) radiation 
emitted by BH-WH pairs entirely comes from the discrete set of modes. 
Using the above equations, the energy carried by <j) of Eq. (Il2| is given by 



E = (<t)\H(t)) = / dujuj (\a UtU \ 2 + K, n | 2 ) + £ a =i,iv HA a 6 a c* + h.c.) . (15) 
Jo 

Because of the complex frequency modes, it is unbound from below. Notice also that the absence of terms such as 
|6 a | 2 is necessary to have at the same time complex frequency eigenmodes and real energies. 



B. Asymptotic behavior and roots k\ 

The material presented below closely follows that of [8.. In fact the lengthy presentation of that work was written 
having in mind its applicability to the present case. The novelties are related to the fact that, for the metrics of 
Eq. ([2|, the supersonic region is bound (from — L to L), and the velocity is subsonic for x —> ±oo. 

Since the velocity v is asymptotically constant for \x/L\ ^> 1, in both asymptotic regions, the solutions of Eq. ^ 
are superpositions of four exponentials e %kxX weighted by constant amplitudes. To characterize a solution, one thus 
needs to know (on one side, say on the left) the four amplitudes associated to the corresponding asymptotic wave 
vectors k(X). These are the roots of 

(\-v(x)k) 2 = k 2 + ^ = n 2 (k), (16) 

evaluated for v(x — >• — oo) = V-. We shall not assume a priori that A is real. Rather we shall look for all ABM. 



Notice that when considering complex frequencies A = u + ir, the roots of Eq. (16) are continuous functions of T. In 
addition, when the scales are well separated, i.e. when k/A <C 1, the relevant values of T will obey T/A <C 1. It is 
therefore appropriate to start the analysis with A = uj real, and then to study how the roots migrate when T increases. 

When uj > 0, since the flow is subsonic for x —> ±oo, there exist two real asymptotic roots: k% > and k% < 
which correspond to a right and a left mover respectively. There also exist a pair of complex conjugated roots, since 



Eq. (16) is real. Thus, on each side of the BH-WH pair there is a growing and a decaying mode. As in [8 we define 



them according to the behavior of the mode when moving away from the BH-WH horizons. 



We are currently trying to demonstrate this assumption. 



Figure 2: Graphical resolution of Eq. (16). On the left, the straight lines represent uj — kv+ in the supersonic regime v+ < —I 
for three different values of uj. The middle one, which is tangent to Q(k), determines the critical frequency cj ma x- For uj < cj ma x, 
the two extra negative roots correspond to right moving modes with respect to the fluid since dkQ > 0. On the right, k u is the 
posi tive root describing the right moving mode, see Eq. (18), the most negative root, k^ that with the smallest norm, see 
Eq. (p22l. 



In preparation for what follows, we study the roots in the supersonic region between the horizons. For uj smaller 
than a critical frequency u; max , whose value is discussed below, the four roots are real. For flows to the left, v < 0, the 
two new real roots correspond to two right movers (with respect to the fluid). Indeed, they live on the u branch of the 
dispersion relation Eq. (16), that with dk^l > 0, see Figj2j When uj increases at fixed \v\ > 1, these roots approach 
to each other. Thus for flows characterized by a maximal velocity u + , there is a frequency u; max above which they no 
longer exist as real roots. It is given by the value of uj where they merge for v = v+. When the two horizons are well 
separated, v+ is given in Eq. (|3|, up to exponentially small terms. As shown in [8], o; max is of the form cj max = Af(D), 
where D is defined in Eq. J2h. For D <C 1, one finds f(D) oc D 3 / 2 . Thus, for a given dispersion scale A, cj max can be 
arbitrarily small. This will be important when considering the appearance of the laser effect in parameter space. 



C. The continuous spectrum 



We now have all the elements to show that the continuous part of the spectrum of ABM is labeled by positive real 
frequencies lj, and that, for a fixed uj, its dimensionality is 2. The general solution of Eq. ^ with uj real can be 
characterized by the 4 amplitudes, which multiply the four exponentials evaluated in the asymptotic left region. When 
imposing that the growing mode is absent on that side, only three independent solutions remain. However, when 
propagating these solutions in the asymptotic region on the right side of the BH-WH pair, the growing mode will be 
generally present on that side. Thus when requiring that it be also absent imposes to take particular combinations, 
and this reduces the dimensionality of ABM to two. 

One can then take appropriate linear combinations to construct the in (out) modes describing the left and right 
movers propagating toward (escaping from) the BH-WH pair. The in right moving solution (/>% m is the combination 



which on the left is asymptotically proportional to e %k ^ x where is the asymptotic real positive root of Eq. (|16). 
Similarly one can identify the in left moving mode </>^ m , and the two out modes (j)^ out and ^ out . 

More precisely, because of the infinite and flat character of the space on either side of the BH-WH pair, the two 
in and the two out modes can be normalized as in a constant velocity flow (by considering a series of broad wave 
packets localized in one asymptotic region and whose spread in frequency progressively vanishes). Considering 



for x — >• — oo, Eq. (13) and Eq. rt9j) imply that it asymptotes to 



^ V ) V duj v^fe) 

A similar expression valid for x — ^ oo gives These two in modes are orthogonal to each other, establishing the 



Kronecker 5 %J in Eq. (13). For uj > 0, these modes have a positive KG norm. For uj < 0, they have a negative norm. 
Therefore negative norm modes can all be described as superpositions of complex conjugated positive norm modes. 

When propagated across the BH-WH geometry, the in modes are scattered by the gradients of v(x). When the 
variation of v is slow, i.e. d x In v <C d x In fc w , the exact solutions are globally well approximated by the WKB solutions 
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of Eq. ([7]). For the u-mode, one finds, see Eq. (17), 

We use the symbol ip™ (ip^) to differentiate the WKB solution from the exact one ^ (</>^). At high frequency, 
uj/k ^> 1, the inequality d x \nv <C d x ln/c^ is satisfied, and u and v modes do not mix. At lower frequency, they do. 
Nevertheless, far away from the BH-WH pair, since v(x) is asymptotically constant, exact solutions decompose into 
superpositions of (p™ and (p^ with constant amplitudes. This applies for both the real and the complex frequency 



modes of Eq. (12). We shall use this fact several times to characterize the properties of the exact solutions. 



Introducing the out modes as in Eq. (|17|), this scattering is described by 

(19) 

Unitarity imposes |T W | 2 + |i? w | 2 = 1 = |T W | 2 + \R U \ 2 , and R^T* + T w i?* = 0. For all values of cj, one thus has an 



elastic scattering, without spontaneous pair creation. For frequencies uj > o; max , Eq. (19) coincides with what is found 
in single horizon scatterings [8 . Instead, for frequencies < uj < u; max , this radically differs from the scattering on a 
single because, in that case, the matrix was 3x3 and mixed </>^, </>^ with the negative frequency ix-mode (4"-lj)*- The 
presence of the second horizon therefore 'removes' these modes. As we shall later see, they shall be 'replaced' by a 
finite and discrete set of complex frequency modes. This is not so surprising since the classical trajectories associated 
with the negative frequency modes are closed (hence the discretization) , and since these trapped modes mix with the 
continuous spectrum through each horizon (hence the imaginary part of the frequency). In fact this is reminiscent, 
but not identical, to quasinormal modes [19] or resonances. The main difference is that quasinormal modes are not 
asymptotically bound. Thus, they should not be used in the mode expansion of Eq. ( [T2| ). 

It should be also mentioned that both in and out modes contain a trapped component in the supersonic region 2 
which plays no role as far as their normalization is concerned since the modes are everywhere regular and the supersonic 
domain is finite. The case where the subsonic domain is also finite should be analyzed separately. If periodic conditions 
are imposed at the edges of the condensate, one obtains discrete frequencies and resonances effects [3 [7]. If instead 
absorptive conditions are used, the frequencies are continuous and the situation is closer to the case we are studying. 



D. The discrete spectrum 

On general grounds we explain why, when considering Eq. ^ in BH-WH flows of Eq. there exists a discrete 
and finite set of complex frequency modes. To this end, we first show that for a generic complex frequency \ = uj + iT, 
there exists no ABM. When r«A, the four roots have not crossed each other with respect to the case where T = 
since the imaginary part of the two complex roots k^j with uj real is proportional to A. Thus, in this regime, one can 
still meaningfully talk about the two 'propagating' roots fc^, fc^, and the growing and decaying roots k^. Then, as in 
the former subsection, when imposing that the growing mode be asymptotically absent on both sides of the BH-WH 
pair, the space of solutions is still two. However when T > 0, the w-root k u acquires a positive imaginary contribution, 



which means that the i^-mode of Eq. (17) diverges for x — » — oo. To get a bound mode, its amplitude should be set 
to 0. For similar reasons, on the right, the incoming v mode diverges for x — » oo. Hence, for a general value of A, the 
set of ABM is empty. The above reasoning applies to flat backgrounds with v constant, and establishes that in that 



case, complex frequency modes should not be considered in Eq. (12). 

We should now explain why when v is supersonic in a finite region, some complex frequency ABM exist. The 
basic reason is the same as that which gives rise to a discrete set of bound modes when considering the Schroedinger 
equation in a potential well (or the propagation of light through a cavity). In the supersonic region, (in the well), 
there exist two additional real roots k u when A = uj real. The classical trajectories associated with them are closed, 
and, as in a Bohr-Sommerfeld treatment, the discrete set of modes is related to the requirement that the mode be 
single valued and bound. The complex character of the frequency A is due to the finite tunneling amplitude across the 



Because of this trapped wave, the in (and out) modes are not asymptotic modes in a strong sense since a wave packet made with <p^ tn 
will have a double spatial support f or t — > — oo: the standard incoming packet coming from x = — oo, and the unusual trapped piece. 
This additional component, see Eq. (I22J, ensures that <$^ %n is orthogonal to the complex frequency modes of Eq. (Ilgj). 
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horizons. Indeed, as we shall later see, the imaginary part of A a is proportional to some quadratic expression in the 
P Bogoliubov coefficients characterizing the scattering through the horizons. Were these coefficients equal to zero, A a 
would have been real. When 'turning on' these coefficients, the frequency A a migrates in the complex plane, and the 
bound modes are continuously deformed. 

These ABM appear in pairs with complex conjugated frequencies. This stems from the hermit icity of H which 
guarantees that there exists an ABM of frequency A* whenever there is one of frequency A a . At this point it should 
be re-emphasized that the existence of these complex frequency ABM is due to the fact that the scalar product is not 
definite positive. Indeed, these modes all have a vanishing norm in virtue of Eq. (11). We can also conclude that the 
discrete set of complex frequencies ABM is finite. In fact there are no closed orbits for uo < — cj max , since the extra 
real roots k u no longer exist and since there is a gap between the eigenfrequencies. 



E. The quantization 

The canonical quantization of the field <p is straightforward since each eigenfrequency sector evolves independently 



from the others. Indeed when decomposing the field as in Eq. (12), with the coefficients a, 6, c promoted operators, 
the equal time commutation x), 7r(t, x')\ = iS(x — x') and the orthonormality conditions Eq. ( [l3| ) and Eq. (14) 
entirely fix their commutation relations. For real frequency modes, the operators o W) j = obey the standard 

commutation relations 

[a w>i , alj] = Sij S(u - J). (20) 
Instead, for complex frequency modes, one gets 

Kcl]=i5 aa ,. (21) 

All the other commutators vanish. 

Because of this disconnection, the ground state of the real frequency modes is stable and in fact subject to no 
evolution. Hence the number of quanta of these modes is constant. The evolution of the states associated with 
the complex frequency modes is also rather simple and described in the Appendix. What remains to be done is to 
determine the properties of the asymptotic fluxes. To this end we need a better understanding of the modes. 



IV. THE PROPERTIES OF THE MODES 

From Eq. ([7]), it is not easy to determine the complex frequencies A a and the properties of the modes. Several routes 
can be envisaged. One can adopt numerical techniques. We are presently modifying [20] the code used in [8j [9] to 
address this problem. One can also bypass the calculation of the eigenmodes and directly compute the propagation of 
coherent states, or the density-density correlation function [21], using the techniques of [22]. This is currently under 
study [23]. One can also envisage to use analytical methods by choosing the flow v(x) as in [16]. This method is also 
currently under study [24]. 

In what follows, we use an approximative treatment which is valid when the two horizons are well separated. Doing 
so, we make contact with the original treatment [1] |F based on wave packets. More importantly, we determine 
algebraic relations which do not rely on the validity of our approximations. In particular, we establish that the real 
frequency modes </>^ are intimately related to the complex frequency modes even though their overlap vanish. 



A. The limit of thin near horizon regions 

To simplify the mode propagation, we assume that the near horizon regions are thin and well separated, i.e. Ln ^> D 
in Eq. In this case, the propagation through the BH-WH geometry resembles very much to that through a cavity. 
Indeed, the following apply. First, the nontrivial propagation across the two thin horizon regions can be described by 
matrices that connect a solution evaluated on one side to that on the other side. Second, the modes can be analyzed 
separately in three regions: in L, the external left region, for — (x + L) ^> D/ Kw] in R, the external right region, 
for (x — L) ^> D l ' kb\ and in the inside region /, for L — \x\ ^> D / n. Within each region, the gradient of v is small. 



Hence, any solution is well approximated by a superposition of WKB waves Eq. (18) with constant amplitudes. 

To further simplify the analysis, we use the fact that the u-v mode mixing coefficients are generally much smaller 
than those mixing the negative frequency modes to the positive u ones [H[9]. Hence, it is a reliable (and consistent) 
approximation to assume that the v modes completely decouple. After having analyzed this case, we shall briefly 
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present the modifications introduced by relaxing this hypothesis. Adopting the hypothesis that the u — v mixing can 
be neglected, for each uj real, one has the following situation. In the left region L, one only has the WKB mode cp™ 



of Eq. (18). Thus, the only solution is 0™' m (up to an overall irrelevant phase we take to be 1). In the inside region 



/, one has three modes: 

c in = a,<pz+ B^H^ir + B^H^ir, (22) 



since in supersonic flows, there exist two extra real roots in Eq. ( |16| ). The superscripts u, (1), (2) characterize the 
coefficients and the WKB modes associated with the three roots shown in Fig. [2] Since we are considerin g a solution 
with idt = oj > 0, the (positive norm) negative frequency modes ip^ appear complex conjugated in Eq. (22). 



In the external R region the solution must be again proportional to the WKB mode <p™. By unitarity, the solution 
must be of the form 0^' m = e %d ^p^. Thus, a full characterization of <^ %n requires to compute the phase 6^ and the 
above three coefficients. At this point, it should be noticed that A = uj is a priori real. However the S matrix, and 
therefore the three coefficients, are holomorphic functions in A. Hence nothing prevents to leave the real axis. In fact 
we shall show that the complex frequencies correspond to poles associated with a coefficient of the S matrix. 

B. An S matrix approach 

1. The S matrix 

The simplest way we found to compute the above coefficients is to follow the approach of [5], up to a certain point. 
In this treatment, a solution of frequency uj > is represented by a two component vector (^,^-u;)- The time 
evolution of a wave packet of such solutions is then considered in the thin horizon limit. Since the frequency content 
of the wave packet plays no role, we do not need to introduce a new notation to differentiate it from an eigenmode. 
In this language, the S'-matrix characterizing a bounce of the trapped mode p^ can be decomposed as 

S = U 4 U 3 U 2 U 1 . (23) 

The first matrix describes the scattering across the WH horizon. In full generality we parameterize it by 

Ul= (a aU z* ^a^) (24) 

Unitarity imposes that \a u \ 2 = \a u \ 2 and |a w | 2 (l — |z w | 2 ) = l. The second matrix describes the free propagation (i.e. 
without backscattering) from the WH to the BH horizon 

In a WKB approximation, the two phases are respectively given by the actions 

J-L 

" k£\x)] , (26) 
In U2, S^l is multiplied by —i since it governs the evolution of y?* w . Its momentum is k™^ = —k^ > 0, where k^ is 



= I dx 



the most negative root of Eq. (16) found in supersonic flows. The ends of integration L w , R u are, respectively, the left 



and right locations of the turning points of the trajectories with — uj < 0. These obey Hamilton's equations [251 ES] 

dx/dt = (dukZ)- 1 , dk u /dt = -d x uj = -k u d x v, (27) 

for negative frequency. In the thin horizon approximation, the turning points hardly differ from — L and L. Unitarity 
brings no conditions on these phases. 

The third matrix describes the scattering across the WH horizon and we write it as 

u "(^kX")- (28 » 
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Unitarity imposes that = 1%^ 2 and ^^(l — |u; a; | 2 ) = l. The fourth matrix describes the return of the negative 
frequency partner towards the WH horizon, whereas the positive frequency mode propagates away in the R region. 
This is described by 



e 



(2) 



(29) 



In the WKB approximation, this backwards movement (hence the +i in the front of S^) is governed by 



1 = J l " dx[-kj?Hx)], 



s& = 



(30) 



(2) 

where the momentum is the least negative ix-root of Eq. ( 16 ). Since the positive frequency mode further propagates 
to the right, there is no meaning to attribute it a phase in U4. In any case this phase would drop from all physical 
quantities. 



The matrix S of Eq. (23) is unitary since its 4 constituents are. Hence IS22I 2 
622 and S21 we shall later use are given by 



\Su 



1 + I S12 1 • The components 



S22 
S21 



-i(S. 



(1) 



i(S. 



(1) _o(2) 



(31) 



Hitherto we followed the method of [5 j. Henceforth, we proceed differently by adding a key element. We require that 
the mode propagated by S be single valued. For real lj, this unequivocally defines (j)™ of Eq. (22). Moreover, when 
looking for complex frequency bound modes, this will give us the modes (p a ,^a we are seeking. 



2. The real frequency modes 
Imposing that the trapped mode of negative frequency is single valued translates into 



i" = s (i ) • < 32 > 



The phase 6 U is that mentioned after Eq. (22). It should not be constrained since the positive frequency component 



keeps propagating to the right. The matricial equation gives 

7 £21 



1 — £22 



Sn 1 — S22 

$22 1 _ ^22 



(33) 



These equations constitute the first important result of this section. They do not rest on the WKB approximation. 
Of course, this approximation can be used to estimate the 4 elements of S. But once these are known, e.g. using a 
numerical treatment, these equations apply. What is needed to get these equations is the neglect of the u — v mixing, 



and no significant frequency mixing in the inside region in order to obtain well-defined amplitudes in Eq. (22). 

Using = S one verifies that the norm of the right hand side of the second equation is unity. This is as it is must 
be, since in the absence of u — v mixing, the ^-component only acquires a phase, here measured with respect to the 
WKB wave. Adopting the convention that WKB modes tp u have a vanishing phase at x = —L, the coefficients of 



Eq. (22) are 



*A W ol^ ( 1 H~ z^j b u ) , 

4 1} = &M + b u ), =6 W . (34) 

These amplitudes are governed by b ui which can a priori be larger or smaller than unity. In particular it diverges if 
622 — >• 1 for some cj, thereby approaching a resonance, see Fig. 3. We now show that the complex frequency modes 
correspond to these resonances. 
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Figure 3: By making use of numerical results borrowed from [20], we have represented the square norm of the amplitude 



solves S22 = 1, see Eq. J36L |££ 2) | 2 behaves as a Lorentzian 



of the negative frequency trapped mode of Eq. (22), as a function of uj real. Near a complex frequency \ a — oo a — iT a which 



see Eqs. (33 34). The dots are the numerical 



values, whereas the continuous line is a fitted sum of Lorentzians. The remarkable agreement establishes that the complex 
frequencies A a can be deduced from the analysis of the real frequency modes. The frequency uj has been expressed in the units 
of Wmax, see Fig. 2, so that there is no resonance above uj/uJm^ = 1. In the present case there are 13 resonances. The narrow 
peaks, T a ~ 0, are due to the fact that the surface gravities are equal, Kb = kw, which leads to z u = Woj in Eq. (41 ). 



3. The pairs of complex frequency modes 



Following the discussion of Sec jIII D| we impose that the amplitude of the incoming u branch be zero, and, as above, 
that the trapped mode of negative frequency is single valued. This gives 



t ).s ; 



which implies 



Pa 



'12, 



s- 



22- 



(35) 



(36) 



Two important lessons are obtained. First, to get the complex frequencies A a with a positive imaginary part, it suffices 
to solve the roots of S22 = 1- Second, as mentioned, these correspond to the poles characterizing the propagating 
modes, see Eq. (33). 



Before computing these frequencies, we explain how to get the decaying modes ip a , the 'partners' of the cp a in 
Eq. (12). Since these bound modes have a negative imaginary frequency, the amplitude of the escaping mode must 



be zero, i.e. 



S 



Using the hermitian conjugated S^, and the unitarity relation S^S = 1, this condition gives 



Pa 
1 



St 







(37) 



(38) 



from which we get (3 a = (6^)12 and 1 = (8^)22- As explained in the Appendix, when expressed in terms of the 
elements of 5, these equations give 



h = -[5 2 i (A*)]*, 1 = [5 22 (A*) 



(39) 



Thus the solutions of 1 = 5| 2 are the complex conjugated of those that solve Eq. (36), thereby establishing the 
partnership between cp a and ip a . In addition, one has f3 a = —(detS)*/3 a . 
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C. The set of complex frequencies A a 



To compute the roots of Eq. ( 36 ) we need to know S22 as a function of A 



+ iT. In what follows we shall relate 
them to the quantities which enter in Eq. (|31|. To this end, we suppose that the 'tunneling' across the horizons is 
small, i.e. the /3-Bogoliubov coefficients associated with each horizon are small. This is true for uj / n sufficiently large, 
see [8] where it was shown that z u and w^, which are related to the Bogoliubov coefficients by z u = j8 w /a w , behave 
as ~ e -*W« (1 _ cj/cjmax) 1 / 4 where cj max is defined in Fig. [2] To proceed we suppose that the norms of w 2 and 
z^Wu are much smaller than 1. In this case, one can expand Eq. (31) in these three products, and in T, since the 
roots of Eq. (36) are real when z u = = 0. 

Indeed, to zeroth order in these quantities, one gets 0^7^ = e~ zn . The norm is trivially constrained by unitarity, 
whereas the phase takes its value from the contributions stemming from the prefactor of the trapped mode at the two 
turning points, see Eq. (18). Taking this into account, 622 = 1 gives 



S. 



(i) 



? (2) 



dx[-k£\x) + k ( J\x)} +ir = 2irn, 



(40) 



with n E Ni. This is the Bohr-Sommerfeld condition applied to the negative frequency mode (In fact, subtracting 

from both and k^ the value of k at the turning points, fc*f , the differences k^ — fcjf , kffl — k 1 ^ have opposite 
sign. Hence, Eq. (40) contains a sum of two positive contributions, as in the Bohr-Sommerfeld condition.) We call 
cj a , a = l,2,...,7V the discrete set of frequencies, which is finite because no solution exists above cj max . 
To first order, for each u a , one gets a complex phase shift 5\ a = 8uo a + %Y a . The imaginary shift is 



The phase in the cosine is 



+ |w Wo | 2 +2|z Wo Wa, | COS^a = ^^(^a)!' 



= S% a + S^l a + arg (z Ua w^ a a UJa /a UJa ) 



and > is the time for the negative frequency partner to make a bounce. It is given by 



d_ 



S. 



(2) 



s 



(1) 



arg a^u 



(41) 



(42) 



(43) 



evaluated for uj = uj a . The first two terms give the classical (Hamilton- Jacobi) time, whereas the last one gives the 
contribution from the scattering coefficients (which is a small correction when LAD 1 / 2 ^> l). 3 The last equality of 
Eq. (41) tells us that T a is linearly related to the norm of the effective /3-Bogoliubov coefficient of the pair, which 

It is also worth noticing that S12 fixes the amplitude /3 a of the leaking mode in E 



obeys^S'22 



|2 



Because | 1 2 = 1 + 1^12| 2 7 Y a defined in Eq. ( |4l| ) is positive, thereby implying that (p a , the solution of Eq. (|35|7is a 
growing mode in time, and an ABM in space. What distinguishes the present case from usual resonances characterized 
by a decay rate (r < 0) is the fact that the norm of the trapped mode y?* w is opposite to that of the leaking wave 
(p™- Even though unitarity in both cases implies a decrease of the norm of the trapped mode, in the present case it 
becomes more negative whereas in the standard case it tends to zero since it has the same sign as that of (p™. As 
a corollary of this, the fact that resonances have an opposite sign of T while satisfying an outgoing condition as in 
Eq. (35) implies that they are not ABM, and therefore not included in the set of modes of Eq. ( [l2| ). This remark 
applies to fermionic fields and implies that the set of ABM for the fermionic dispersive field considered in pQ and 
propagating in the BH-WH metric of Eq. ( [2|, i s restricted to the continous set of Eq. (12), i.e. positive real frequency 
modes elastically scattered as given in Eq. ([19]). To complete these remarks, one should notice that ip a: our decaying 
modes, are also ABM because they obey the incoming condition Eq. ( [38] ) . 

Our treatment is similar to the interesting analysis presented in [11] . In that work, the general solution is constructed 
in terms of WKB waves. Then the Bohr-Sommerfeld and the outgoing conditions are separately imposed in the small 
tunneling approximation, thereby fixing both the real part of the frequency uo a and its imaginary part T a . We followed 
another logic which leads to the same result, namely, the requirement that the mode be an ABM gives Eq. (35), which 
in turn gives the complex equation 622 = 1 that encodes both conditions. We note that S22 — 1 does not require that 
the tunneling amplitudes be small to be well defined. We also note that our quantization scheme applies to the cases 
studied in [HI [12] and allows to remove the 'formal trick' used in the second paper. 



3 Using the techniques of 26 28 , we verified that it gives a small correction to the classical actions. We also obtained a refined expression 
for Eq. GuT) which goes beyond the WKB approximation and which has been validated by the numerical analysis of |20| . 
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D. The density of unstable modes 

To further characterize the instability, it is of interest to inquire about the density of unstable modes, about the 
end of the set, its beginning, and about the most unstable mode. 

The number of unstable modes will be either large, or small, depending on the value of L/uj mdiX . When L/uj max ^> 1, 
the number is large because, for u < cj max , the gap between two neighboring modes roughly given by D/L. The end 



of the set is controlled by u; max as explained after Eq. (40). It does not significantly contribute to the instability when 
ft/^max <^ 1, since and w u tend to zero as e _7ru; / K (1 — c^/^max) 1 ^ 4 for uj — » cj max . Thus the growth rate r a — » 
asw^ cj max . 

The beginning of the set is governed by the first Bohr-Sommerfeld modes. Their frequency is of the order of D/L. 
A proper evaluation of T a is harder since, generically, the coefficients a u and j u both diverge as ~ cj -1 / 2 . We thus 



conjecture that the most unstable mode is the first mode, or one of the first ones, because the cosine of Eq. (41) could 
in some cases lower the value of Ti below that of the next ones. We can also characterize the migration of the roots 
as L increases. The (total) variation of Eq. ( |40| ) with respect to L and Eq. ( |43| ) tells us that duj a /dL > 0. From this 
we conjecture that new unstable modes (obtained by increasing L) appear with uo a = 0. This is corroborated by the 
first runs of [20 . The determination of T a for uj a — >• is difficult, and further study is needed to establish if T a follows 
the abrupt behavior present in Fig. 1 of [T4] . 

E. Taking into account the u — v mixing 

A nonvanishing u — v mixing hardly modifies the above results. However the algebra becomes more complicated 
because the two nontrivial matrices describing the propagation across the WH and the BH horizon (Ui and Us) are 
now the 3x3. Let us only sketch this enlarged case. 

When considering the in mode </>^' m , in the left region L, one has </>^' m = ip™ + R^cp^^ see Eq. (19). In the inside 
region /, one has 

4> u J in = Al V l + Al & + + B^i^lT, (44) 

and in the external R region, the mode asymptotes to T^cp^. As before, the coefficients are determined by the 
matching conditions at the horizons. Expressing m in terms of three coefficients associated with the three WKB 
waves, ((/?^, <P-u))i using the Bogoliubov matrices [8] Swh, Sbh which give, for the WH and the BH respectively, 
the coefficients after the scattering given the incident ones, one has 



Swh \ K , (45) 





and 

= S BH I I . (46) 

In the second equation, the exponentials arise from the propagation from the WH to the BH horizon, and from our 
choice that the phase of WKB waves vanishes at the WH horizon. The six relations fix the six coefficients. T w and 




Ru characterize the asymptotic scattering, see Eq. (19), whereas the four others determine <^ %n in the inside region. 



To get the complex frequencies A a , since the weight of the incoming v mode is already set to 0, one simply needs 



to replace the 1 in the right hand side of Eq. (45) by a 0, as done in Eq. (35). 



V. PHYSICAL PREDICTIONS 

A. Classical settings 

From an abstract point of view, any (bounded) initial condition, i.e. the data of the field amplitude and its 



derivative with respect to time, can be translated into the coefficients a^,6 a ,c a of Eq. (12). This is guaranteed by 



the completeness of the mode basis. Yet, it is instructive to study more closely how initial conditions describing a 
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wave packet initially moving towards the WH-BH pair translate into these coefficients. This will allow us to relate 
our mode basis to the wave packet analysis of [l] |5] . Moreover it is a warming up exercise for the determination of 
the fluxes in quantum settings. For simplicity, we work in the regime in which the number of unstable modes is large, 
i.e. Lcj max ^> 1, and with a wave packet of mean frequency obeying uj <C cj max and Luo ^> 1, e.g. with uj = y^uJm&x/L. 

We consider a unit norm wave packet which is initially in the left region and propagating to the right. At early 
times, before it approaches the WH horizon, it can thus be expressed in terms of the WKB modes of Eq. (18) as 



poo 
JO 



As an example, one can consider the following Fourier components 

p -(uj-ou) 2 /2o~ 



e iujt e iS™(x ,-L) 



(47) 



(48) 



The first factor fixes the mean frequency Co and the spread a taken to satisfy cj/uo <C 1 and aL ^> 1. The last two 
exponentials ensure that at t = to the incoming wave packet is centered around some xq <C —L. S%(x,x f ) is the 
classical action of the u-mode evaluated from x to i', see the phase of <p™ in Eq. (18). 
When decomposing ^ as in Eq. (12), the coefficients are given by the overlaps 



k = h) ± o, c a = h) = o. 

Hence, at all times, the wave packet can be expressed as 

POO 

C(M) = / <^ (/., I0 H e""V: in (^) + c-c) + S a (b a e- iA °VaW 



c.c. 



(49) 



(50) 



The first two coefficients in Eq. (49) are easily computed and interpreted. They fix the real frequency contribution 



of the wave packet. The last two are very interesting and encode the instability. We first notice that, because of 
Eq. (14), they are given by the overlap with the partner wave. Hence the coefficients c a identically vanish. Indeed 
the overlap of our wave packet with the growing modes cp a vanishes since these, by construction, have no incoming 
branch. On the contrary b ai the amplitudes of the growing modes, do not vanish since the decaying modes ijj a contain 
an incoming branch, see Eq. (38). It is thus through the nondiagonal character of Eq. (fl4|) that the instability enters 



in the game. In this respect, our analysis differs from Ref. [16]. We do not understand the rules adopted in that work. 

To compute b a we use the fact that at time to the wave packet is localized around xq <C —L. Thus we need the 
behavior of ip a in this region. Using the results of Sec. |IV C and Eq. (18) extended to complex frequencies, for 
T a T* < 1, one gets 



iPa(x) = (3 a ipl a (x)=(3 a ( P Z a (x)xexp[-r a t UJa (x,-L)], 



(51) 



where t UJa (x, — L) = d 0J S u {x 1 —L) > is the time taken by a ix-mode of frequency uj a to travel from x to — L. 4 Inserting 
this result in (^ a |^), using the fact that the wave packet has a narrow spread in x — xq given by l/(ad u k u )^ the 
overlap is approximatly 



0* 



7T 1 /4(2a) 1 /2 



-i(uj— u a )t e iS™ 



(x ,-L) 



x e 



-T a t Ua (x Q -L) 



(52) 



The meaning of this result is clear. The spatial decay of ijj a governed by T a has the role to delay the growth of 
the amplitude of the wave packet until it reaches the WH-BH pair. The phase ensures the spatial coherence of the 
propagation. That is, the sum of a in Eq. (47) will give constructive interferences along the classical trajectory 



emerging from xq at to. Let us also mention that, because of Eq. (14), the normalization of b a is arbitrary, but the 
product b a ip a is well defined and physically meaningful. 

From this analysis, we see that at early times the wave packet ^ behaves as described in [TJ[5]. It propagates freely 
without any growth until it reaches the horizon of the WH. Then one piece is reflected and becomes a v mode with 



4 It is interesting to notice how ip a and <j) a acquire a vanishing norm. Because the leaking wave has an amplitude S21 oc T a ' and decreases 
in x in r« 1 , its positive contribution to the norm is independent of T a and cancels out that negative of the trapped mode. 
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amplitude The other piece enters in the inside region. A later times, a component stays trapped and bounces 
back and forth while its amplitude increases in as e rt . Every time it bounces there is leakage of a u wave packet at 
the BH horizon, and a v mode at the WH one. 

This back and forth semiclassical movement goes on until the most unstable mode, that with the largest T a , 



progressively dominates and therefore progressively destroys the coherence of the successive emissions. In Sec jlVD 
we saw that the most unstable mode is likely to be that characterized by the smallest real frequency, n = 1 in Eq. (40). 
This implies that at late time, our wave packet will be completely governed by the corresponding wave (fi(x) (unless 
of course b\ = 0) 



e 



ri(t-to) 



Rele-^hip^x)}. (53) 



In the late time limit, its behavior differs from the above semiclassical one. Indeed, in the inside region, one essentially 
has a standing wave whose amplitude exponentially grows (it would have been one if the tunneling amplitudes were 



zero). Outside, for x ^> L, using <pi(x) ~ /3i(p™ (x) and the same approximation as in Eq. (51), one has a modulated 
oscillatory pattern given by 



-» e 



riKt-to)-^ x Re[e- <W1 * 61)81^(2;)]. (54) 



It moves with a speed equal to ~ 1 in the dispersionless regime, when cji/A <C 1. The energy flow is now given by a 
sine squared of period 2i\juj\ rather than being composed of localized packets separated by the bounce time T b . 



B. Quantum settings 

1. The initial state 

Because of the instability, there is no clear definition of what the vacuum state should be. Indeed, as can be seen 



from Eq. (15), the energy is unbound from below. Therefore, to identify the physically relevant states, one should 
inquire what would be the state, or better the subset of states, which would obtain when the BH-WH pair is formed 
at some time to. If this formation is adiabatic, the initial state would be close to a vacuum state at that time. That 
is, the expectation values of the square of the various field amplitudes would be close to their values in minimal 
uncertainty states, with no squeezing, i.e. no anisotropy in (j)^ — tt^ plane, where ir^ is the conjugated momentum 
of (fiu. Because of the orthogonality of the eigenmodes, this adiabatic state would be, and stay, a tensor product of 
states associated with each mode separately. 

There is no difficulty to apply these considerations to the real (positive) frequency oscillators which are described 
by standard destruction (creation) operators a™, (a^,a^). The adiabaticity guarantees that one obtains a state 



close to the ground state annihilated by the destruction operators. Because of the elastic character of Eq. (19), for 
all real frequency oscillators, one gets stationary vacuum expectation values with no sign of instability. 

There is no difficulty either for the complex frequency oscillators described by b a and c a . Indeed, as shown in the 
Appendix, one can define two destruction operators <i a+ and d a - , and use them to define the state as that annihilated 
by them at to. In this state, we get the following nonstationary vacuum expectation values: 



<Mt)St(0> 



$a>a 



e -iu a (t-t') e r a (t+t'-2to) 



<<v(*)cj(f)> = ^ e --.(*-*') e -r.(*+*'-2to) } 

{k'{t)c\{t')) = i*|£ e -».(*-0. (55) 

As expected, the expression in bb^ (cc^) leads to an exponentially growing (decaying) contribution, whereas the 
cross term is constant at equal time. This behavior is really peculiar to unstable systems. Even though the metric is 
stationary, there is no normalizable state in which the expectation values of the 6, c operators are constant. Stationary 
states do exist though, but they all have an infinite norm, see the Appendix. 
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Using Eq. (12), and putting to = for simplicity, the two-point function in this vacuum state is 

Ef=iRe(e-^(*-*VaWC(^) 
^-*'V (zK(z')) 

- i(A «*- A: *'Va(x)^*(x') 

*(^*- A »*'V a (x)^*(ai')) . (56) 

(57) 



We notice that the last term, the second contribution of complex frequency modes is real, as a classical term (a 
stochastic noise) would be. We shall return to this point below. We also notice that the second term is purely 
imaginary and, when evaluated at the same point x = x' , it is confined inside the horizons since ip a vanishes for 
x <C — L whereas ijj a does it for x ^> L. Thus it will give no asymptotic contribution to local observables. 



2. The asymptotic fluxes 

Our aim is to characterize the asymptotic particle content encoded in the growing modes ip a . To this end it is useful 
to introduce a particle detector localized far away from the BH-WH pair. We take it to be sitting at x ^> L, in the 
R region on the right of the BH horizon. We assume that it oscillates with a constant frequency cuo > 0, and that its 
coupling to (j) is switched on at t = — oo, and switched off suddenly at t = T ^> to in order to see how the response 
function is affected by the laser effect a finite time after the formation of the BH-WH pair at to = 0. 

When the detector is initially in its ground state, to second order in the coupling g with the field, the probability 
to find it excited at time T is given by [27] 



P e (T) = g 2 f dt' f dte-^-t'Ufax)^^)), 

J — oo J —oo 

= g 2 Z a \/3 a \ 2 K»| 2 | / dte-^°- x ^\ 

Jt UJa {L,x) 

(w - w„) 2 + (r„) 



9 ^7- " , v» n Ua (T,x), (58) 



where 

n Ua (T,x) = |/3 a | 2 e ^ a (T-t^(L, x )) ^ _ exp{ _[r a + i{uQ _ Wa)](T _ ^j} |2_ (59) 



To get the second line of Eq. (58), we used cp a (x) ~ f3 a (p™ a (x), the fact that the BH-WH pair is formed at t = 0, 
and that it takes a time t 0Ja (L^xJioi the mode y?^ to reach the detector at x. To get the third line, we used the 



inequality T a T^ a < 1 as in Eq. (51). The meaning of the various factors appearing in P e is the following. The sum 
over a means that all unstable modes contribute, but the Lorentz functions restrict the significant contributions to 
frequencies uj a near co>o, that of the detector. The prefactor \(p^ a (x)\ 2 depends on the norm of the corresponding mode 
evaluated at the detector location, as in the usual case. The function n 0Ja (T, x) acts as the number of particles of 
frequency uj a received by the detector at time T, and at a distance x — L ^> k^ 1 from the BH horizon. It depends on 
the number initially emitted (= \P a \ 2 ) multiplied by the exponential governed by T — t UJa1 the lapse of time since the 
onset of the BH-WH pair minus the time needed to reach the detector at x. 

From the response function of a localized detector, it is clear that one cannot distinguish between the noise due to 
quanta of the real frequency modes </>^ and that carried by the growing modes ip ai because both modes asymptote 
to the WKB waves cp™ which are asymptotically complete. In this respect it is particularly interesting to compute 
the de-excitation probability Pd which governs the (spontaneous + induced) decay of the detector. It is obtained 
by replacing ujo by — ujo in the first line of Eq. ( [58] ) [27] . In that case, one finds that the spontaneous decay only 
comes from the (f)™ whereas the induced part only comes from the cp a . The induced part equals that of P e since the 
asymptotic contribution of the (p a to Eq. (|57|) is real. We are not aware of other circumstances where orthogonal modes 



16 



with different eigenfrequency (here <fi™ and cp a ) are combined in this way in the spontaneous + induced de-excitation 
probability P^, or equivalently, contribute in this way to the commutator and the anticommutator of the field, i.e. 
with the if a only contributing to the latter. (For damped modes, the commutator and the anticommutator are related 
differently, see e.g. Appendices A and B in [30.) The lesson we can draw is the following: even though the modes 
(j)™ are orthogonal to the growing modes (p a , their respective contribution to the asymptotic particle content cannot 
be distinguished by external devices coupled to the field. 

It is also interesting to compute the asymptotic outgoing energy flux (T uu (t,x)), where T uu = [(d t — d x )(j)] 2 . At 
large distances in the R region, using cp a ~ Pa ( P\ a and Eq. (57), the renormalized value of the flux is 



(T uu (t,x)) = E a \(dt-d x )e-^^ a (x)\' 

= £ a |/? a | 2 |(d t -d x )e- iA «VA>)l 



(60) 



It only depends on the discrete set of complex frequency modes. Yet, because of the imaginary part of A a defines a 
width = T a , the spectrum of real frequencies uo is continuous, as can be see in Eq. (58). In fact the observables can 



either be written in terms of a discrete sum over complex frequencies, or as a continuous integral of a sum of Lorentz 
functions centered at uo a an d of w idth T a . However this second writing is only approximative and requires that the 
inequality T a T^ <C 1 of Sec. IV C be satisfied to provide a reliable approximation. 



3. The correlation pattern 



As noticed in pQ, because of the bounces of the trapped modes, the asymptotic fluxes possess non-trivial correlations 
on the same side of the horizon, and not across the horizon as in the case of Hawking radiation without dispersion [2TJ 
[27| [29] , or in a dispersive medium [25 . These new correlations are easily described in the wave packet language of 
that reference, or that of Sec. V A When using frequency eigenmodes, they can be recovered through constructive 
interferences, as in Sec. IV F. of [9 . Indeed, when the complex frequencies modes form a dense set so that the 
dispersion of the waves can be neglected, the sum over a in Eq. (57) constructively interferes at equal time for two 
different positions x and x' separated by a propagation time 



d UJ S u (x,x f ) 



dx d/jj k^j (^c) 



(61) 



where is the bounce time of Eq. (43). This is because the differences uj a — uj a +n are equal n x 2n/T^ since the i 



are solutions of Eq. (40). With more precision, the conditions for having these multimodes interferences are, on one 
hand, cj max /ft ^> 1 so that dispersion hardly affects the modes and, on the other hand, kL ^> 1, so as to have many 
modes for uo below the Hawking temperature ~ n. (For higher frequencies |/3 a | 2 , which governs the intensity of the 
correlations, is exponentially damped.) 

However, since dispersive effects grow and since the most unstable mode progressively dominates the two-point 
function of Eq. (57), at sufficiently large time the above multimode coherence will be destroyed and replaced by the 



single mode coherence of the most unstable one. This is unlike what is obtained when dealing with a single BH or 
WH horizon because in that case [29] the pattern is stationary, and all frequencies steadily contribute (significantly 
for uu < k). In the present case, at late time, if <pi is the most unstable one, the correlation pattern is given by 



(4>(t,x)4>(t',x')) = 



ri(t+t') 



x Re 



(e-^C-'V^KOr') 



(62) 



The asymptotic pattern, for x ^> L, is obtained using (fi(x) ~ f3\e r i t ^i( L ^ x ) ( pU^^ x y j^. - 1S ver y s i m ii ar to that of 



Eq. (54) found by sending a classical wave packet, see Appendice C of [9 for a discussion of the correspondence 
between statistical correlations encoded in the two-point function and deterministic correlations encoded in the mean 
value when dealing with a wave packet described by a coherent state. 

So far we worked under the assumption that the u — v mixing coefficients are negligible. When taking them into 
account, one obtains a richer pattern which is determined by the complex frequency modes solutions of Sec. |IVE| 



4- The small supersonic region limit 

When L of Eq. ([2| (or D = \v+\/c — 1) decreases, the number of solutions of Eq. (40) diminishes. Therefore, in 
the narrow supersonic limit nL — » 0, there is a threshold value for nL given D, below which there is no solution. In 
that case, there are no unstable mode, and no laser effect. In fact no flux is emitted, and this even though the surface 
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gravity of the BH (and that of the WH) is not zero. The reason is that there is no room for the negative frequency 
modes (j)-^ to exist. In agreement with the absence of radiation, the entanglement entropy of the BH [[31] would 
vanish, because it accounts for the number of entangled modes across the horizon and thus of opposite frequency, see 
[32] for the effects of dispersion on the entanglement entropy. 



5. Comparison with former works 



It is instructive to compare our expressions to those obtained in pQ and in [5]. Our expressions differ from theirs 
because the discrete character of the set of complex frequency modes was ignored in these works. As a result, a 
continuous spectrum was obtained. Yet this spectrum possesses rapid superimposed oscillations stemming from the 
interferences that are present in S2i{oo) of Eq. (31). A priori one might think that they could coincide with our 



(2) 

frequencies uj a = ReA a . However, as noticed in [5 , their value are insensitive the phase governed by Si , whereas it 



In the limit of Sec. IV C 
-cj. Moreover (p { 

where T b , 



plays a crucial role in Eq. (40). We found no regime in which the two sets could approximatly agree. Therefore, as 
far as the fine properties of the spectrum are concerned, the predictions of [TJ [5] are not trustworthy. 

Nevertheless, when the density of complex frequency modes is high, and when ignoring these fine properties, the 
average properties derived using [HIS] coincide with ours. Indeed when considering the mean flux in frequency intervals 
Au ^> 1/L, the rapid oscillations found in (TJ |5] are averaged out. As a result the mean agrees with that over the 
contributions of complex frequency modes. This can be explicitly verified by comparing the norm of our discrete 
modes (p a ,i>a with the continuous norm of the negative frequency modes used in [U [5] . 
relevant contribution to the overlap (ifj a \(p a ) comes from the negative frequency mode (p 

given by a sum of (normalized) WKB waves and times a prefactor = 

time given in Eq. ( [43] ) . This is just what is needed for approximating the discrete sum in Eq. (57) by a continuous 
integral J duo with a measure equal to 1 . 

In conclusion, when the number of bound modes becomes small, the difference between our description and the 
continuous one increases. This difference is maximal when there are no solution of Eq. (40). In this case, no radiation 
is emitted, something which cannot be derived by the continuous approach of [TJ[5]. 



the 
are 



and ip, 
is the bounce 



VI. CONDITIONS FOR HAVING A LASER EFFECT 



Having understood the black hole laser effect, it is worth identify in more general terms the conditions under which 
a laser effect would develop. We define a 'laser effect' by the fact that a free field possesses complex frequency ABM 



while being governed by a quadratic hermitian Hamiltonian, as in Eq. (10), and a conserved scalar product, as in 
Eq. ([9]). The field thus obeys an equation which is stationary, homogeneous, and second order in time. Let us note 
that this type of instability is often referred to as a dynamical instability [6j[T6j|34], a denomination which indicates 
that quantum mechanics is not needed to describe/obtain it. Let us also note that we do not consider the case where 
the frequency of the ABM is purely imaginary. Such dynamical instability seems to belong to another class than that 
we are considering, see the end of this subsection for more discussion on this. 

Using semiclassical concepts, the conditions for obtaining complex frequency ABM are the following : 

• For a finite range of the real part of frequency cj, WKB solutions with both signs of norm should exist, or 
equivalently, positive norm WKB solutions should exist for both sign of uj. This is a rather strong condition 
which requires that the external field (gravitational or electric) must be strong enough for this level crossing to 



take place, i.e. for the general solution be a superposition as in Eq. (22). 



• These WKB solutions of opposite norm must mix when considering the exact solutions of the mode equation. 
In other words there should be connected by a nonzero tunneling amplitude. This is a very weak condition as 
different WKB branches are generally connected to each other. In our case it means that z u and appearing 



in Eq. (31) should not vanish. 



• One of these WKB solutions must be trapped so that the associated wave packets will bounce back and forth. 
This is also a rather strong condition. 



The depth of the potential trapping these modes should be deep enough so that at least one pair of bound 



modes can exist, see Eq. (40). This condition is rather mild once the first three are satisfied. 



When these conditions are met for a sufficiently wide domain of frequency cj, they are sufficient to get a laser effect, 
and they apply both when the external region is finite [6] or infinite. Being based on semiclassical concepts, stricto 
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sensu, they cannot be considered as necessary. But we are not here after mathematical rigor, rather we wish to 
identify the relevant conditions in physically interesting situations. 

In this respect, it should be noticed that when only the first two conditions are satisfied, one obtains a vacuum 
instability [27 , also called a superradiance in the context of rotating bodies [33j [34]. Hence, whenever there is a 
vacuum instability, one can engender a laser effect by introducing a reflecting condition, as was done in [T3l [14] , or by 
modifying the potential, so that the last two conditions are also satisfied. It should be clear that when the laser effect 
takes place, it replaces the vacuum instability rather than occurs together with it. Indeed, as proven in Sec. |IV C[ 
the frequency of the trapped modes are generically complex. The possibility of having a trapped mode (subjected 
to a vacuum instability prior introducing the reflecting condition) with a real frequency is of measure zero, as two 
conditions must be simultaneously satisfied. To give an example of the replacement of a vacuum instability by a laser 
effect, let us consider the archetypal case of pair production in a static electric field studied by Heisenberg [35] and 
Schwinger [36 . In that case, one obtains a laser effect by replacing the Coulomb potential Ao = Ex by Aq = E\x\ 
which traps particles of charge q for qE < 0, for frequencies uo < — m where m is their mass. We hope to return to 
such pedagogical examples in the near future. 

In conclusion, we make several remarks. Even though the complex frequency ABM are orthogonal to the real 



frequency modes, as is it guaranteed by Eq. (11 ), the asymptotic quanta associated with these modes are not of a new 
type but are, as we saw, superpositions of the standard ones associated with real frequency modes. If laser instabilities 
can be studied in classical terms, the quantum aspects are not washed out. For instance, when considering a charged 
field, the charge received as infinity is still quantized, albeit its mean value is described by a complex frequency ABM. 
Moreover, in all cases, when the instability ceases, the number of emitted quanta is a well-defined observable governed 



by standard destruction/creation operators as those appearing in Eq. (A13). From this it appears that dynamical 
instabilities governed by an ABM with a purely imaginary frequency, as e.g. the Gregory-Laflamme instability [37] . 
belong to another class since this asymptotic decomposition in terms of quanta does not seem available. Whether it 
could nevertheless makes sense to quantize such instability is a moot point. 

VII. CONCLUSIONS 

We showed that the black hole laser effect should be described in terms of a finite and discrete set of complex 
frequency modes which asymptotically vanish. We also showed that these modes are orthogonal to the continuous set 
of real frequency modes which are only elastically scattered, and which therefore play no role in the laser effect. In 
Sec. |IV[ using the simplifying assumption that the near horizon regions are thin, we determined the set of complex 
frequencies and the properties of the modes using the approach that combines and generalizes [5] and [11] . 

We described how an initial wave packet is amplified as it propagates in the BH-WH geometry. When the density 
of complex frequency modes is high we recovered the picture of [1] at early times. Instead, at late time, or when the 
density is low, the successive emissions of distinct wave packets associated with the bouncing trajectories are replaced 
by an oscillating flux governed by the most unstable mode. 

We then computed in quantum settings how the growth of the complex frequency modes determine the asymptotic 
fluxes when the initial state at the formation of the BH-WH pair is vacuum. Because of the width associated with 
the instability, the spectral properties of the fluxes are continuous albeit they arise from a discrete set of modes. The 
properties we obtained significantly differ from those found in [TJ[5]. We also briefly described the properties of the 
correlation pattern at early times when the number of complex frequency modes is large, and at late time when only 
the most unstable mode contributes. When the supersonic region between the two horizons is too small so that there 



is no solution to Eq. (40), we concluded that there is no instability, that no flux is emitted, and that the entanglement 



entropy vanishes. Finally, in Sec. [VI] we gave the general WKB conditions under which a laser effect would obtain 
starting from the standard concepts that govern a vacuum instability in Quantum Field Theory. 

This work poses several questions which deserve further study. In Bose condensates, the backreaction due to the 
instability and the suppression of the instability could be computed using the Gross-Pitaevskii equation. Fermionic 
fields jU [12] should be further studied to reveal the roles of the N pairs of complex frequency modes, which are not 
ABM, but which indicate that the naive vacuum will decay in the lowest energy state vacuum plus N asymptotic 
quanta by spontaneously emptying the N Dirac holes which are trapped inside the horizons. 
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Appendix A: Upside down harmonic oscillators 

a. Real upside down oscillator 

We review the quantization of upside down oscillators. To begin with, we start with a single real upside down 
harmonic oscillator. Its Hamiltonian is 

tf = ^(p 2 -rV), (Ai) 

when written in terms of position q and conjugated momentum p, obeying the standard equal time commutator (ETC) 
[q,p] = i- Introducing the 'null' combinations 

1 (p + rg), c=-L(p-r g ), (A2) 



V2T V2T 

one gets 



H= ^(bc + cb). (A3) 



One verifies that they obey the ETC [6, c] = i. The ordering of b and c in H follows from that of Eq. (Al). The 
equations of motions are 

b=(-i) [6, H] = T 6, c=(-i) [c, H] = -T c, ( A4) 

thereby establishing that b (c) is the growing (decaying) mode b = bo e rt , (c = Co e~ vt ). 

It is now relevant to look for stationary states. In the ^-representation (c = — idb), the stationary Schroedinger 
equation H^e = E^e reads 

-ibd b * E (b) = (^- l -)* E (b). (A5) 

Solutions exist for all real values of E, and the general solution is 

*E(b) = A E 0(b)(b) lE/T - 1/2 + B E 0(-b)(-by E/r - 1/2 . (A6) 
Since the spectrum is continuum, one should adopt a Dirac delta normalization (^e'\^e) — ${E — E f ). This gives 

i^r+is B i 2 = 2^f. ( A? ) 

Imposing that the solution be even in q (p, 6, or c) imposes Ae = Be- The important lesson one should retain is that 
there is no square integrable stationary states. Therefore, in all physically acceptable states (i.e. square integrable) 
the expectation values of q 2 + p 2 will exponentially grow ~ e 2Ft at late time. 



b. Complex oscillators 

More relevant for us is the complex upside down harmonic oscillator. It can be described by the complex variables 
q = qi + iq2, p = Pi + ip2 where qi,pi and ^2,P2 are hermitian and obey the standard ETC given above. We then 
introduce the complex b and c variables 

1 (p + Tq), c=^=(p-Tq), (A8) 



/AT V4r 
which are normalized so that they obey the ETC 

[6,c f ] =i. (A9) 
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We then look for the (hermitian) Hamiltonian which gives the following equations 

b = (-i)[b,H] = -iXb, c=(-i)[c,H] = -i\*c, 

where A = to + iT. It is given by 

H = -i\Jb + i\*tfc. 



(A10) 



(All) 



It is instructive to reexpress this system is term of a couple of destruction and creation operators d_, dl, d + , d+ which 
are given, at a given time to, by 



b = -j=(d+ - idi), c = -j=(-id+ + dl). 



(A12) 



They obey the standard commutation relations [d^,dj] = 5ij, and Eq. (All) reads 



H = u{d\_d+ -did-) +T(d_d + + d 1 f dL), 



(A13) 



In the first term one recovers the standard form of an Hamiltonian is the presence of stationary modes with opposite 
frequencies [9 . The second term induces a squeezing of the state of the d_ , d + oscillators which grows linearly with 
time. 

To set initial conditions, and to be able to read the result of the instability in terms of quanta, it is appropriate to 
use this decomposition of H and to work in the 'interacting' picture where the operators 6, c only evolve according to 
Ho, and where the squeezing operator acts on the state of the field. Indeed, in this picture the states can be expressed 
at any time as a superposition of states with a definite occupation numbers n_ and n + . 

c. The hermitian conjugated 



To make contact with the treatment of Sec. IV B, it is appropriate to express the time evolution of linear operators 
in a S- matrix language. We introduce the operator 



A[Z{t)\ = x{t) b + y{t) c = (b co) , 



(A14) 



where bo and cq are the operators 6, c evaluated at to- Using Eq. (A10), the time evolution of Z = (x,y) is 



where the 2x2 matrix H is diag(A, A*). By definition, the S matrix brings Z from to to to + 1. It is given by 

-i\t q 



e 



-iX*t 



(A15) 



(A16) 



To define we need to refer to the matrix encoding the ETC of b and c, see Eq. (A9): 



Q 



[6,6t] M] 
[c,ftt] [c,ct] 



i 

-i or 



This matrix defines a scalar product in the space of the vectors Z through 

[A(z'),A(zy] = (z\z') = {x* y*)Q^ yl )- 

The hermitian conjugate defined by (SZ\Z r ) = (Z\S^ Z') obeys 

q S ] = S * T Q. 



(A17) 



(A18) 



(A19) 
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In the present case, using Eq. (A16), we find 



5t(A) = 







JX*t 



= [S(A*)]*- 



(A20) 



In the case where S would have been nondiagonal, we should also transpose the matrix: 

5 f (A) = [S(A*)]* T . 

We see that the expression of the components of 5* in terms of those of S i s unu sual because the norm of Z, instead 
of being \x\ 2 + \y\ 2 , it is of the form lm(xy*). However, when A is real, Eq. (A20) gives the usual relation. Hence, for 



(A21) 



complex A, it can be viewed as an analytical continuation of the real frequency case. 

The hermicity of the Hamiltonian guarantees that the S matrix is unitary in the sense that the scalar product {Z'\Z\ 
or the ETC relations Eq. (A9), are conserved. One verifies that Eq. (A20) implies the usual unitarity condition 



^ (A) 5(A) 



(A22) 



for A real and complex. 

The above analysis applies to the modes ijj a when considering the field operator restricted to the a sector, i.e. 
A(Z a ) = xb a (p a + yca^a, because ip a , ifj a have complex conjugated frequencies and because b a , c a obey Eq. (21). 
Hence, the hermitian conjugated matrix appearing in Eq. (38) should be understood as in Eq. (A21). 
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